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ABSTRACT 



The Palomar Transient Factory (PTF) Orion project is one of the experiments within 
the broader PTF survey, a systematic automated exploration of the sky for optical tran- 
sients. Taking advantage of the wide (3.5° x 2.3°) field of view available using the PTF 
camera installed at the Palomar 48-inch telescope, 40 nights were dedicated in 2009 
December to 2010 January to perform continuous high-cadence differential photometry 
on a single field containing the young (7-10 Myr) 25 Ori association. Little is known 
empirically about the formation of planets at these young ages, and the primary mo- 
tivation for the project is to search for planets around young stars in this region. The 
unique data set also provides for much ancillary science. In this first paper we de- 
scribe the survey and the data reduction pipeline, and present some initial results from 
an inspection of the most clearly varying stars relating to two of the ancillary science 
objectives: detection of eclipsing binaries and young stellar objects. We find 82 new 
eclipsing binary systems, 9 of which are good candidate 25 Ori or Orion OB la associ- 
ation members. Of these, two are potential young W UMa type systems. We report 
on the possible low-mass (M-dwarf primary) eclipsing systems in the sample, which 
include six of the candidate young systems. Forty-five of the binary systems are close 
(mainly contact) systems, and one of these shows an orbital period among the shortest 
known for W UMa binaries, at 0.2156509 ±0.0000071 days, with flat-bottomed primary 
eclipses, and a derived distance that appears consistent with membership in the general 
Orion association. One of the candidate young systems presents an unusual light curve, 
perhaps representing a semi-detached binary system with an inflated low-mass primary 
or a star with a warped disk, and may represent an additional young Orion member. 
Finally, we identify 14 probable new classical T-Tauri stars in our data, along with 
one previously known (CVSO 35) and one previously reported as a candidate weak-line 
T-Tauri star (SDSS J052700. 12-^010136.8). 

Subject headings: Binaries: close - binaries: eclipsing - open clusters and associations: 
individual (25 Ori) - planets and satellites: detection - stars: pre-main sequence - 
techniques: photometric 



1. Introduction 

The Palomar Transient Factory (PTF) is a survey built around a wide-field mosaic camera 
installed on the Palomar 48-inch Samuel Oschin telescope dedicated to surveying the sky to find 
photometric transient and variable objects with variability timescales of minutes to years. The 



^Einstein Fellow 
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camera consists of a mosaic of two rows of six 2kx4k-CCDs (of which one is non- functional) , with 
1" pixels, givin g a total held siz e of ~ 3.5° x 2.3°. With a 5a detection limit of i? ~ 21.0 mag for 
60s exposures (jLaw et al.ll20ld ). the main survey is fully automated, and is designed to allow for 
very rapid followup of transient events using the Palomar 60-inch telescope, as well as a further 
global n etwork of te l escop es for additional characterization. A tech nical summary of the system is 
given by lLaw et al.l (120091 ). and an overvi ew of the science goals by lRau et al.l (120091 ) . A summary 
of the first year's performance is given by lLaw et al.l (|2010l ). 



Among the programs within the survey is the PTF Orion project. The primary goal of the 
Orion project is to search for transiting exoplanets around stars that are 5-10 Myr old. Little 
is known empirically about exoplanets during the first few millions of years of their lives - the 
timescales on which they are expected to form and migrate - and the aim is to tr y and fill this ga p 



in order to provide constraints on theories of planet formation and evolutio n (e.g., 
Typical young circumstellar disk lifetimes are on the order of 5 — 10 Myr (jHillenbrand 



Armitagd 20091 '). 



20081) ■ and 



by trying to catch stars at or just after the point of disk dissipation, we can hope to find planets 
at the point where they may first become observable without their photometric signatures being 
swamped by that of the disk, and where disk-planet interaction has ceased. 

In addition to the planet search, the observations provide a unique data set to study a variety of 
stellar astrophysics: eclipsing binary systems enable tests of star formation and evolution models; 
stellar activity and rotation at these young ages can be characterized; and previously unknown 
pre-main-sequence (PMS) stars can be identified and characterized wit hin the region. (A si milar 
young trans i t and eclipse search, the Monitor project, is described by lAigrain et al.l ()2007l ) and 
Miller et al.l l{20o4 ).) 



The 25 Ori association, with an estimated age of around 7-10 Myr (|Briceho et al.l 120071 ). pro- 
vides an ideal target to achieve the goals of the PTF Orion survey, being the most populous known 
sample in its age range within 500 pc, the right age, and with a spatial extent we ll matched to 



the PTF fie ld of view. The field has already been well studied and characterized by iBricefio et al 



(|2005l . 120071 ). In the winter of 2009/2010, we obtained time-series data centered on the 25 Ori as- 
sociation, taking continuous i?-band observations and obtaining ~ 110,000 individual high-cadence 
light curves. 

The unique capacity of eclipsing binary systems for yielding accurate measurements of fun- 
damental stellar parameters, such as mass, radius, and effective-temperature ratio, makes them 
useful for constraining stellar models. This is especially true o f young PMS syste ms, where 
evolutionary models are not as well constrain ed observationally (jMathieu et al.l 120071 ) . At pri- 



mary masses below 1.5 Mq, iHebb et al.l (|2010l ) list in addition to thei r newly reported system 



MML 53, o nly six known PMS echpsing binaries (RXJ 0529. 4+0041A, 



Irwin et al 



V1174 Ori. IStassun et al. 2004; 2MASS J05352184-0546085, Stassun et al 



Covino et al 



20071: Pa r 1802, iCargile et al.l 120081 : IStassun et al.l 120081 : and ASAS J052821+0338.5 



good . 



20061 . l2007l : JW 380 



2004 



Stempels et al.l |2008| ) , with a further four candidates in the Orion Nebula cluster more recently 
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reported by iMorales-Calderon et al.l (|201ll ). Evolutionary models for stars at the late end of the 
stellar mass range (^MO) are also currently not well constrained by empirical measurements, owing 
to their intrinsic faintness, and the consequent dif ficulty in observing them. To date, fewer than 
20 such systems are known (see iKraus et al.l l201ll . and references therein) . It is therefore of par- 
ticular interest to identify binary systems which are young Orion members and/or which lie in the 
low-mass regime in our data set. 

Observations of contact binary systems, commonly known as W UMa systems, are also of in - 
terest both because they are potentially useful distance indicators (lRucinskill2004l : lEker et al.ll2009l) 



and b ecause their formation and evolution is cur rently not well understood (see, e.g., iBilir et al 
20051 : lEker et al.l l2006l : iPribulla &: Rucinskil 120061 ) . It has been debated whether they may form 



as the result of the merger of previously detached binar y systerns thr o ugh a ngular momentum 
loss due to tidal interaction and magnetic braking (e.g., IStepiehl Il995l . l2006l ). or as a result of 
initially detached binary orbits decaying by a combinatio n of the Kozai mechanism and tidal 



friction due to the presence of a third m ore distant body (jTokovinin et al.l l2006l : iRucinski et al 



20051 : [Roxburgh 



196 



a. 



2007 : lEggleton &: Kisseleva-Eggletonll2006l ). or instead directly as contact systems (e.g.. 



Bilir et al 



The 7-10 Myr age of the 25 Ori region, for exam ple, is too short for 



the expected Gyr timescale needed for in -spiral via angular momentum loss (jStgpienlll995l . 



Bilir et al 



2005 



2006 



Gazeas fc Stgpiehl 120081 ) . and indeed, appears somewhat short even for merger 



by the Kozai /tidal- friction mechanis m, which operates on timescales of tens of Myr or more, 
(jEggleton &: Kisseleva-EggletonI |2006| ) . Finding young contact binaries within our sample is po- 
tentially useful for constraining these formation mechanisms. 

Single PMS T-Tauri stars (TTSs) are also useful for understanding stellar formation and early 
evolution. At the age of the 25 Ori association, they are relatively unobscured and c an be observed 



in th e optical regime. The advent of wide-field surveys like the CIDA-QUEST survey ([Briceho et al 



20051 ) allows for broad, deep searches that are sensitive to the more dispersed and more mature 
low-mass populations, where previous studies have been more tightly focused on the brighter OB 
stars in relatively compact regions. The PTF survey is similarl y sensitive, an d identificati on of new 
young TTSs in our data complements the sample reported by lBriceho et al.l (|2005l . 120071 ). 



In this first paper we present a description of the survey and the data reduction procedure, 
and some of the initial results that relate to these two stellar astrophysical science objectives of the 
program: binary systems and new candidate PMS stars. Section [2] describes the survey; section 
[3] describes the data-reduction pipeline developed for the PTF Orion project; section H] gives an 
overview of the data obtained and an assessment of the current precision levels; section [5] discusses 
the eclipsing binaries found, including the candidate young binary systems, the low-mass binaries, 
and some of the contact and near-contact binaries of interest (with some overlap between these 
groups); and section [6] discusses the new classical T-Tauri stars (CTTSs). A brief summary is given 
in section [3 
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2. The PTF Orion Survey 

2.1. Field Selection 

Our target field around 25 Ori was chosen to balance a number of requi rements. Interstella r 
reddening is minimal along the line of sight to the cluster (mean Ay = 0.29; iBricefio et al.l 120071 ): 
the source count density is optimal, being high but without overcrowding (~ 4.4 sources per square 
arcminute at our detection thresho ld); and the region i s already relatively well characterized, with 
an estimated age of ~ 7-10 Myr (jBricefio et al.l 120071). correspond ing to the time at which the 
disks around young stars have mostly dissipated (jHillenbrandl l2008l ) . This allows us to probe the 
photometric variability of young stars without it being completely dominated by variations due to 
accretion or circumstellar extinction. The majority of TTSs in the 25 Ori region are weak- lined T 
Tauri stars (WTTSs), indicating that the disks have indeed dissipated (jBriceho et al.l 120071 ). 



Within these constrai nts, the field was aligned to maximize the number of known TTS covered 



(Bricefio et al 



2005 



20071 : C. Briceho 2009, private communication) whilst allowing for the non- 



functioning chip, and to place 25 Ori itself on a gap between chips, minimizing charge bleeding from 
what would otherwise be a heavily saturated source. Figure [T] shows the positioning of the field 
with respect to 25 Ori and the brighter stars in the region. The field is centered at approximately 
a = 05^'26™46?, 6 = +01?50f50'.' (J2000), covering an area 3.5° in right ascension and 2.3° in 
declination. One of the chips in the array is non-functional (the north-east of the central four, chip 
number 4) and is not shown. 

The figure confirms that the amount of interstellar reddening is apparently unifo rmly quite low 
in the 25 Ori field, and we adopt the derived mean extinction value for the cluster from lBriceho et al] 
(|2007l . Ay = 0.29 mag) as a singl e representative value. (In comparison, typical values from the 
dust maps of lSchlegel et al.l (119981 ) along lines of sight to the sources in the field range from « 0.3 
to 0.6 mag, and mostly < 1.0 mag. As estimates of the full line-of-sight extinction through to the 



edge o f the Galaxy, these suggest an upper limit.) Following the extinction laws of ICardelli et al 
(jl989l ). and assuming Ry = 3.1, the adopted value of Ay corresponds to extinctions in R, J, H, 
and Ks of 0.22 mag, 0.082 mag, 0.055 mag, and 0.033 mag, respectively. 



2.2. PTF Observations 

Data were taken for the Orion project on the majority of the clear nights between 2009 De- 
cember 1 and 2010 January 15. Weather allowing, the chosen primary field was observed in the R 
filter as near continuously as possible, whenever it was higher in the sky than an airmass of 2.0. 
Exposures were 30s long, with a cadence varying generally between 70-90s, including readout time 
and depending on the performance of the telescope control system and guiding control loop (see 
below) . 

Since the Orion observing program could not easily be incorporated into the normal PTF 
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Fig. 1. — : Location of the field with respect to 25 Ori and the Orion belt stars, overlaid on grayscale 
contours representing mean 2MASS J ~ Kg colors, which trace reddening. Bin size is 0.1° x 0.1°. 
Filled circles mark bright stars {V < 6.0 mag), with size varying linearly with visual magnitude; 
the brightest of these are labeled, along with 25 Ori which is placed in a chip gap in the PTF 
field of view. The footprint for each chip in the field is marked and the chip number identified; 
th e dead chip (NE of th e central four) is omitted. The circle shows the 1° working radius used 
bv lBriceno et al.l (|2007l l to define the 25 Ori group. 




160 170 180 190 200 210 220 

JD - 2455000 

Fig. 2. — : Number of exposures taken per night for the PTF Orion project from December 2009 - January 
2010. Solid black regions indicate the fraction of exposures which remained unflagged after 
processing by the data reduction pipeline, for chip (similar fractions are seen for the other 
chips) . 

robotic scheduling, separate scripts for telescope and camera control code were adapted from the 
standard PTF robot code for the purposes of the Orion observations. Each night at the beginning 
of Orion field observations, normal PTF operations were manually interrupted to run the Orion ob- 
serving program. When the field center set once again below airmass 2.0, control was automatically 
returned to the standard PTF robot scheduler and normal PTF operations resumed. 

Of the 40 nights dedicated to the Orion project, 14 yielded usable data, the majority of the 
remainder being either of marginal quality due to cloud or subject to telescope closure because of 
poor weather. (One night also could not be processed due to an insufficiency of non-Orion exposures 
to create a flat field for that night, and can in principle be recovered.) A total of 4,486 exposures 
were taken during the run. An average across the chips of ~ 3, 460 exposures ran successfully 
through the pipeline after rejection of those which could immediately be discarded at the outset by 
visual inspection, and an average of « 2, 400 of those exposures passed the data reduction pipeline's 
automated data quality assessment without flagging. There is some variation in these numbers from 
chip to chip since each is processed completely independently. Figure [2] shows a histogram of the 
number of exposures taken on each night, and the fraction of those which passed unflagged for chip 
0, which is used as a representative example of all the chips. 

In order to minimize the effects of flat-fielding errors for high-precision photometry, it is desir- 
able to minimize image shift on the detector between frames. Since the Palomar 48" has no native 
guiding facility, we developed software to guide on the science images as they came from the tele- 
scope, providing feedback to the telescope via offset pointing commands. Offsets were referenced 
to a single good quality reference frame taken at the beginning of the run, and telescope pointing 
updates were applied during readout time generally every second frame. On initial pointing of the 
telescope, offsets of typically ~ 70-90 pixels relative to the reference frame were normally corrected 
within two exposures. With the guiding switched on, the pointing toward the center of the field was 
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stabilized with a root-mean-square (rms) image shift of ~ 1.0 and ~ 0.5 pixels in right ascension 
and declination respectively, compared to unguided linear tracking drifts of ~ 0.3 pixel min~^. 



2.3. Ancillary Data 



In addition to the PTF data obtained, we also drew on data from several other sources. The 



photographic US Naval Observatory USNO-Bl.O source catalog (|Monet et al.l 120031 ) provided a 
reference data set for initial photometric zero-point correction of the -R-band PTF Orion data 
(section 13. A further catalog of ~ 40,000 reference stars from the Centro de Investigaciones 
de Astronomia-Quasar Equatorial Survey Team (CIDA-QUEST), provided by C. Briceho (2009, 
private communication), allowed for a final more accurate correction to the i?-band zero point after 
differential photometry had b een performed (section l3.2p . This catalog had in turn been calibrated 
against the Landolt system ([Landoltl Il992l ) with R magnitudes ranging from 14 < < 19 and 
measurement errors varying with increasing magnitude from ~ 0.015 to 0.1 mag. The overall 
error for the brighter stars (down to R ~ 17) includ ing systematics is ~ 0.02-0.03 mag. We also 
made use of the PPMXL catalog (jRoeser et al.ll2010l ) to obtain proper motions for our candidate 
binary systems. The PPMXL catalog provides a determination of positions and proper motions 
by combining USNO-Bl.O and Two Micron All-Sky Survey (2MASS) astrometry, and aims to be 
complete to y ~ 20 mag. 



The 2MASS point-source catalog (jSkrutskie et al.l 120061 ) provides all-sky coverage of point 
sources in the near-infrared J (1.25/xm), H (1.65/im), and Kg (2.16^m) bandpasses, with magnitude 
limits of 15.8, 15.1, and 14.3 mag respectively. PTF Orion source detections were matched against 
this catalog to provide color information on each source in combination with the PTF Orion R 
magnitude. The 2MASS catalog includes a 'quality' flag for each measurement in each band, which 
can take various values. In particular, designations from A' to 'F' represent photometric quality in 
decreasing order (in the case of an 'F', error estimates could not be determined). In certain cases 
throughout this paper we use these flags to restrict analysis to sources with 2MASS data above a 
given quality level, as stated in the text. 

For certain target objects, anticipating that these young objects might have discernible mid- 
i nfrared (MIR) exce sses due to circumstellar disks, we also examined the Spitz er Space Telescope 
(jWerner et al.ll2004l ) Heritage Archive for both Infrared Array Camera (I RAC, iFazio et al.ll2004l ) 
and Multiband Imaging Photometer for Spitzer (MIPS. iRieke et al.ll2004l ) data, discussed below. 
We found some IRAC data (at 3.6, 4.5, 5.8, and 8/im), but MIPS data (at 24/im) were more common 
in our target region. We found MIPS and/or IRAC data for 15 objects, discussed in sections 15.2.51 
andU All 15 of the objects were found in at least one of 9 AORKEYtll|: 10987776, 26066432, and/or 



^AORKEYs are the unique eight-digit identifier for the Astronomical Observation Request (AOR), which can be 
used to retrieve these data from the Spitzer Heritage Archive. 
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26067456 for IRAC, and 10987008, 17052160, 26067712, 26067968, 26069760, and/or 26070016 for 
MIPS. The images for each object were examined for potential source confusion, and none was 
found. 

For IRAC, we started from the automatically-produced mosaics from the Archive (S18.7 or 
greater), and performed aperture photometry using the IDii routine aper.prjf] on these mosaics at 
the designated positions. We used an aperture radius of 3 native pixels (1.2 arcsec pixel"^), a sky 
annulus of radius 3-7 native pixels, and multiplicative aperture corrections of 1.124, 1.127, 1.143, 
and 1.234 for the four channels, respectively, as prescribed by the IRAC Instrument Handbook, 
available on the Spitzer Science Center (SSC) Web site|l| The zero-points we used to convert between 
flux density units (Jy) and magnitudes were, respectively, 280.9, 179.7, 115.0, and 64.13 Jy. 

For MIPS, we initially obtained the automatically-produced mosaics from the Archive (S18.12 
for MIPS), and on a case- by-case basis, reprocessed the mosaics if necessary from the basic cal- 
ibrated data (BCD) level, to reduce the influence of instrumental artifacts. We also conducted 
aperture photometry using IDL's aper.pro on the mosaics (either ours or the Archive's), using an 
aperture of 7", an annulus of 7"-13", and an aperture correction of 2.05, from the MIPS Instrument 
Handbook, again available from the SSC website. The MIPS-24 zero point we used was 7.14 Jy. 

In the ideal case, to determine whether or not there is infrared excess due to a circumstellar 
disk, we need an estimate of a spectral type so that we can fit a model to the spectral energy 
distribution (SED) and estimate reddening (Ay), and then estimate the degree of infrared excess 
due to a disk. In the absence of an accurate estimate of spectral type, the degree of IR excess is 
not well defined. We can instead attempt to estimate the degree of infrared excess by assuming 
that these sources are all truly stars (as opposed to background objects), and then comparing the 
longer Spitzer bands to a shorter band. Where possible, even lacking a spectral type estimate, we 
can use [3.6] — [24] to reasonably reliably indicate disk excess (where the bracket notation denotes 
the magnitudes at the Spitzer bands), since this range is within the Rayleigh- Jeans tail even for 
the latest spectral types, giving them colors of zero. Reddening also does not strongly affect this 
color unless Ay is extreme, in which case the extinction will be high enough to make a source 
detection unlikely in our PTF data. MIPS-24 data are more widely available in our field than 
the IRAC bands, however, and 3.6/um data do not always exist. In such cases, the 2MASS Kg 
band is the most commonly available alternative. Kg — [24] therefore provides an alternative excess 
estimate, although it is le ss reliable since late -type stars (such as many of those in our sample) are 



not colorless at Kg — [24] (jOautier et al.ll2007l ). and the effects of reddening are stronger at Kg than 



at 3.6/im. For such cases, additional followup observations will be needed to refine the estimate of 



^|http: //www, ittvis . coni/idl| 



^Found in the IDL Astronomy User's Library, [http : // idlastro . gsf c . nasa . gov/ 1 see also iLandsmarJ ( 19931 ) and 



StetsonI (|l987l ) 



' jhttp: //ssc . spitzer . caltech. edu/ 1 
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the excess. 



3. The PTF Orion Data Reduction Pipeline 

The initial PTF data reduction pipeline in place at the time the data were taken was geared 
primarily towards absolute photometry and extensive database operations for handling the vast 
quantities of data coming from the broader survey. We developed a separate dedicated relative- 
photometry pipeline, which took over from the standard pipeline at the point where image process- 
ing was complete. Here, we describe the two broad steps in the data reduction procedure: obtaining 
the initial raw photometry, beginning with the standard PTF image processing and the first steps 
of the Orion pipeline; and the differential photometry process used to obtain the precisions needed 
for the PTF Orion survey. 



3.1. Raw Photometry 

The main PTF pipeline image processing consisted of standard bias subtractions, flatfielding, 
and astrometric solutions. Each chip was reduced separately and in parallel. Floating biases 
were measured as the median of the overscan regions of the chips and subtracted along with a 
nightly combined superbias frame; flatfields were created by combining all (non-Orion-project) 
images take n with the same filter du ring a night, all the bright pixels being first masked out using 
SExtractor ([Bertin &: ArnoutsI Il996l ). to leave a p ure sky fiatfi eld. Astrometric information was 
added to the FITS image headers using SCAMP (IBertinl 120061) to solv e for the astrometry. The 
general PTF pipeline is discussed in more detail in lGrillmair et al.l (|20ld ) and R. Laher et al. (2011, 
in preparation). 



The Orion pipeline ran separately and took as input the image data from the standard pipeline 
after the astrometric solutions had been added. Written in IDL, the main steps in the PTF Orion 
pipeline are as follows: 



1. Pixel masking: Masks are initially applied to reject any badly behaved pi xels on the detecto r 
using an algorithm loosely based on the IRAfI ' ccdmask' procedure (jTodvl Il986l . Il993l ). 
Masks are created from images made by dividing a 70s light-emitting diode (LEe[§) fiatfield 
by a 35 s LED fiatfield; three independent such divided frames were obtained for each chip. 
Any pixels with outlier fiuxes beyond 4c7 in at least two of the three frames, or beyond 
3a in all three of the frames were fiagged as bad. This approach helps catch excessively 
variable pixels in addition to non-linear pixels, while still rejecting cosmic-ray hits. The 



^http: //iraf .noao . eduT] 

"^See lLaw et al] f\20oi ) 
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flagging procedure was then repeated several times after boxcar smoothing of the original 
image along the readout direction with a selection of bin sizes from 2 to 20 pixels (this 
removes columns where individual pixels are not statistically bad, but are statistically bad 
when taken as an aggregate). Pixels lying in small gaps between bad pixels were then also 
iteratively flagged, with the aim of completely blocking out large regions of bad pixels while 
minimizing encroachment into good pixel regions. Any source detection with a flagged pixel 
within its photometry aperture is rejected and does not pass through the pipeline. 

2. Raw photometry: Source detection and aperture photomet ry are performed using the ID L 



Astronomy User's Library implementation of DAOPHOT (iLandsmanl Il993l : IStetsonI 119871 ). 
with minor modiflcations to return more significant figures in the output files and to handle 
the large number of sources in each CCD frame. Owing to the sometimes large offsets of the 
field on the detector from image to image, source detection is performed separately on every 
exposure obtained, with a Aa detection threshold. A 4-pixel aperture radius provides a rea- 
sonable compromise between maximizing the fiux enclosed for bright sources and minimizing 
sky-background noise for faint sources, while also being large enough to prevent significant 
photometric errors from imperfect centroiding of the apertures. This compares to median 
seeing of 2.5" at full- width-half-maximum (FWHM; where 1 pixel ~ 1"). Inner and outer 
radii of 10 and 20 pixels, respectively, are chosen for the annulus used for background estima- 
tion around the aperture. The default limits for 'sharpness' and 'roundness' of the detected 
sources provide a first cut to reject cosmic rays and background galaxy detections. 

3. Initial zero point correction: An initial zero-point correction is applied by matching detec- 
tions agamst the USNO-Bl.O catalog and oonect.ng tor the outlier-res.stanfl mean difference 

between raw and USNO-B magnitudes for each image from each chip. This places all pho- 
tometry on an absolute scale, within the accuracy of the USNO-B catalog (~ 0.3 mag). Sloan 
Digital Sky Survey (SDSS) photometry data would provide a more accurate correction, but 
were not available for the 25 Ori field. 

4. Creation of master source list: After filtering out images with poor seeing or high background, 
a master source list for each chip is created by stepping through each of the 100 images with 
the highest numbers of source counts, in decreasing order of source counts. For each image, 
if a source is detected whose sky position cannot be matched with one already found (within 
a 2" radius, the matching radius used throughout the pipeline), it is added to the master list. 
On every 20th image, sources with only one detection are pruned from the table, as are any 
sources with fewer than 10 detections at the end of the process. This ensures that cosmic 
rays, asteroids, and spurious detections are rejected from the table. Since we are primarily 
concerned with persistent stellar objects with relatively low-level variability, we accepted the 
consequent loss of short transient events that may appear from below the detection limit in 



^Throughout this text, where an outlier-resistant mean or standard deviation is referred to, these are found using 
the 'resistant_mean' and 'robust_sigma' functions, respectively, from the IDL Astronomy User's Library. 
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return for the significant savings in data storage space in the final output data tables. This 
does not affect any of the sources that are suitable for planet transit searches. 



3.2. Differential Photometry 

After the initial raw photometry and preparation of the master source lists, a single data table 
for each chip is prepared with entries for each source in the lists at each epoch of observation. A 
second pruning of the table is performed to remove sources if they show very few total detections 
(fewer than 100) and those detections are sparsely distributed in time (the fraction of successful 
detections within any 50 consecutive epochs is less than 0.5). This helps remove sources which 
are due to badly behaved pixels or on the extreme edge of the detection limit, while retaining 
genuine transients that are only detected for a limited length of time. The source list from each 
chip is then split into four equal-area quadrants, each treated separately. This provides a first-order 
accommodation for variations in photometry across the wide field of view which may be caused 
by thin cloud, variations in point spread function (PSF) d ue to the wide-field opt ics, effects of 
differential airmass variation, and atmospheric refraction (see Everett Sz HowelllboOll ). On average 



~ 2600 sources were detected per chip quadrant for the 25 Ori field. 

In a given quadrant, an outlier-resistant mean magnitude is then calculated for each source 
across all epochs; this creates a reference magnitude list against which to perform the differential 
photometry. For each observation, the difference is calculated between every source's measured 
magnitude and its corresponding reference magnitude. Within each quadrant, the mode of the 
differences is then calculated and subtracted from the measured magnitudes to provide the differ- 
entially corrected magnitudes. The mode is calculated by fitting a Gaussian curve to the peak of 
the histogram of magnitude differences in the quadrant, and the peak of this Gaussian is taken 
to represent the best-likelihood estimator of the true fine-scale correction for the zero point of the 
frame. All stars in a given quadrant are thus used in establishing the differential correction applied 
for that quadrant. 

Using the modal magnitude difference correction in this manner has the advantage of being 
very robust to outliers and variable stars. It also requires no selection of an ensemble set of stable 
reference stars, provided that there is a sufficient number of stars in the field and that a sufficient 
number of those stars are relatively stable. Even in the event that only a small fraction of stars 
are stable, provided that the total number of stars is large, the mode still represents a reasonable 
estimate of the differential correction since there is no reason to expect any correlation in intrinsic 
variability between the stars included in the histogram. Additionally, the modal difference draws 
information from the full ensemble of stars, providing a natural 'weighting' of the differences: 
stars that have higher photometric uncertainties (i.e., lower fluxes) or that are more variable will 
contribute a more spread-out Gaussian to the total distribution, adding more to the wings of the 
histogram, though still peaking around the true magnitude difference; those with lower uncertainties 
(high fluxes) will contribute a tighter Gaussian to the total and contribute more to the peak of the 
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histogram. 

To further improve the statistics, the differential photometry is then repeated after identifying 
and omitting from the histogram sources which stih show an rms variation greater than 4(7 above 
their mean measurement error, or an rms variation greater than 2% (or more precisely, an rms 
greater than the quadrature sum of these two values). This entire procedure is iterated until the 
total number of sources rejected in this way changes by less than 3%. 

As a final step, once differential photometry is complete, the real apparent magnitude is set 
by comparing the median magnitudes of each differentially corrected time series against the list 
of stable reference stars from the CIDA-QUEST survey, which in turn are calibrated against the 
Landolt system (see section 12. 3p . The median of the differences is used to provide a final single 
zero-point correction for each entire chip. 

After processing, the data are automatically checked and any questionable sources, nights, or 
individual data points are fiagged according to a number of criteria. Flags are added for failed 
DAOPHOT photometry or differential photometry; multiple source detections, USNO-B reference 
sources, or master-list sources, within the FWHM; high image background, very low zero-point, or 
high variation in background across the image; exceptionally large FWHM; low number of sources 
across the image; large variations with time in the measured right ascension/declination for an 
individual source location; saturation; high sigma or skew in individual background annuli for 
sources; proximity to predicted optical ghosts; and indication of crowding or contamination within 
the photometry aperture (assessed by finding outlier values for the ratio of fiux inside the aperture 
to that in a thin annulus around the aperture compared to other sources in the field). 

4. Data 

The differential photometry pipeline yielded a total of 116,285 individual light curves (including 
some resulting from false multiple detections on saturated stars). For an initial investigation of 
the clearly variable stars, a subsample of all the light curves with an outlier-resistant rms variation 
greater than three times the outlier-resistant mean measurement-error for each light curve was 
selected. This cut alone yielded ^ 2060 sources. These light curves and their corresponding 
images were visually inspected to select only those which clearly represented intrinsic photometric 
variation, rather than crowded sources or image artifacts which had escaped fiagging, sources on 
unmasked bad pixels, or sources whose variability was clearly due to known systematics (see section 
14. ip . This initial inspection yielded a total of ~ 530 sources, the difference in number from the first 
cut being primarily due to the ~ 4mmag noise floor at the bright end which was not accounted 
for, and to biasing toward light curves that were not too sparsely sampled. (For comparison, 
selecting only sources with > 1000 of a possible ~ 2400 unflagged detections, and rms greater than 
three times the quadrature sum of the mean measurement-error and the noise floor, yields ~ 700 
light curves). Among the light curves selected from the initial inspection, we identifled visually 82 
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eclipsing binaries and 16 likely classical T-Tauri star (CTTS) light curves. The results presented 
in this paper are based on these last two groups of light curves. 



4.1. Precision and Sensitivity 

Figure [3] shows a plot of the standard deviation of unflagged data points versus R magnitude for 
all of the ~ 100, 000 light curves obtained, on all 11 chips, over the full length of the observations. 
The results are compared with the theoretical errors including photon, background, and readout 
noise. Some chips are better behaved than others, but the noise floor at around 4mmag is evident. 
Approximately 8% of all light curves with more than 10 detections have an rms < 1%. The banding 
in the theoretical errors comes from the different chips, and is likely due to small differences in 
quantum efficiency. Saturation begins to set in around R ~ 13 mag, where the upturn in the 
rms appears. The faint-end limit varies depending on the chip, but the turn-over in the number 
of source detections with increasing faintness begins to set in around R ~ 19.5 mag, with few 
detections beyond 20.5 mag. 

In addition to intrinsically variable stars, some of the scatter of points above the noise floor at 
brighter magnitudes in figure [3] is due to systematic trends that appear in small regions of certain 
chips. After performing differential photometry, it is not unusual to discover common systematic 
trends in some of the resulting light curves that may be a function of position, color, and/or airmass. 
An advantage of the 'modal differencing' technique employed in the differential photometry (section 
13. 2p is that, provided the fraction of sources in the ensemble exhibiting a given trend is relatively 
small, the position of the peak of the difference histogram should be insensitive to the trend. As a 
result, sources that are not directly affected by the trend will also not be compromised as a result 
of the differential photometry corrections. 

The presence of trends appears to be fairly minimal in the PTF Orion differential photome- 
try, with the exception of one effect which is attributed to known fogging of the detector during 
the period of the Orion observations, due to slow deposition of oil contam inants from the dry-air 



system used at the time to prevent condensation on the CCD window (see lLaw et al.ll2010l ). This 
contributes a slow long-term trend in certain regions of the detector over the length of the obser- 
vations, with the exception of the last night, just before which the CCD window had been cleaned. 
After the cleaning the affected light curves returned to their initial flux level. The trend causes 
variation of up to ~ 0.2 mag in the worst cases, but is generally easy to recognize, and appears to 
be confined to relatively small regions on the detector; only a small fraction of the light curves are 
significantly affected at the level of variability with which we are concerned here. Chip 4 suffered 
the most from the fogging effect, which leads to a slight overdense band of sources from that chip 
around 0.02-0.03 mag in figure [3j At this point light curves which clearly show the trend are re- 
jected from our current analysis for simplicity (four probable eclipsing binaries are affected in this 
way). A number of sources on chip 4 alone also show random noise on the few-percent level on the 
first three nights and the last night of the run, which may be detector noise related. This should 
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Fig. 3. — : Standard deviation (black points) vs. median magnitude, compared to theoretical photometric 
errors (cyan), for all unflagged data from all light curves obtained with more than 10 good data 
points. There is a slight overdense band around 0.02-0.03 mag attributable to chip 4, where a 
known systematic trend affects certain localized regions. See section [4. II 
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be borne in mind with regard to results from that chip. 

There are some suggestions of other smah trends, possibly color related, but no detrending of 
the data has been attempted at this stage, and none of these effects are accounted for in figure [3j 
Nonetheless, the figure shows that the majority of the data are well behaved down to the 4mmag 
noise floor without any detrending. 

4.2. Absolute Photometric Accuracy 

The absolute accuracy of the magnitudes is dependent on the second zero-point correction 
against the CIDA-QUEST reference stars (see section [23]) . The outlier resistant standard deviation 
in the magnitude difference between the CIDA-QUEST sources and the least variable 50% of the 
matched PTE Orion sources ranges from 0.025-0.044 mag depending on the chip. Combining this 
with the overall 0.02-0.03 mag error of the CIDA-QUEST list (including systematics) , we therefore 
conservatively estimate the zero points of the individual light curves to be accurate to w3%-5% for 
chips 0-4 and 6-10. 

The area covered by chips 5 and 11 in the Orion field had no coverage in the reference star set, 
and therefore no direct secondary zero-point correction could be determined. Instead, a correction 
was made based on the weighted mean of the corrections from the other chips, excluding only chip 
4, owing to more significant fogging-trend systematics on that chip which give rise to a distinctly 
bimodal, non-Gaussian distribution in the zero-point histogram for that chip. The mean of the 
secondary zero-point corrections for the remaining chips was 0.159 mag, with a standard deviation 
of 0.044 mag. Adding this whole-chip zero-point uncertainty in quadrature with the typical uncer- 
tainty for the individual light curves discussed above, and a 0.03 mag error for the CIDA-QUEST 
list, yields an absolute accuracy estimate of 0.069 mag for the light curves on chips 5 and 11. 

Marginal color-dependent trends are seen in the magnitude residuals against the CIDA-QUEST 
references (as a function of R—I), but they are below the level of the random scatter in the residuals, 
implying that errors due to mismatch between the PTE and CIDA-QUEST filters are negligible 
at our accuracy level. Since the error estimates for the zero points are based on the scatter in the 
residuals about the mean, rather than the error in the mean itself, any systematic errors in zero 
point should be included in the accuracy estimate. 

5. Eclipsing Binary Systems 

Among the most easily identified of the light curves are the eclipsing binaries. In this section 
we present those found in this initial investigation. Section 15.11 presents an overview of the data; 
an analysis is presented in sections 15.21 where we identify candidate young binaries, and 15.31 where 
we identify some additional low- mass systems. In section [53] we summarize and discuss the specific 
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noteworthy sources. 



5.1. Eclipsing Binary Results 

We identified a total of 82 clear eclipsing binaries in the PTF Orion data, for which period- 
folded light curves are presented in figures SHHl These were identified by visual inspection of 
our short-list of strongly varying light curves, most eclipsing binary light-curve types displaying 
distinctive characteristic shapes. Light curves which were sinusoidal in appearance were rejected 
due to the ambiguity in interpretation of their nature. Am ong the binar i es, th ree were found to 



have been previously identified as candidate variable stars bv lKraus et al.l (|2007l ) in the MOTESS- 
GNAT survey (sources 0-9653, 5-12287, and 5-12446). In those cases where a period could not 
be determined, either because too few eclipses were detected or because the data did not show a 
regular period, the data covering only the best captured eclipse are instead plotted against time to 
show the eclipse shape. The complete light curves for the five such targets are shown in addition 
in figures M - flOl 

Table [1] lists the basic properties of the eclipsing binary systems: the assigned PTF Orion 
survey ID, N-n, consisting of the chip number, N, on which the source was detected, followed by a 
running sequence number, n; the mean measured J2000 coordinates over the whol e observing run 



(accu rate to ~ 1"); the corresponding 2MASS ID where a match was obtained (jSkrutskie et al 



20061 ): the 2MASS J, H, and Kg magnitudes; classification as 'close' (C) or 'detached' (D) (see 
below); the period, P, in days; Tq, the heliocentric Julian date for the epoch of the primary eclipse; 
the estimated distance to the system; the median measured R magnitude; and AR, the approximate 
peak-to-peak magnitude range of the measured light curve. Specific sources discussed in the text 
are broken out at the beginning of the table. The newly measured properties in the table are 
obtained as described below. 



Classification: The 82 systems are categorized into one of two broad classifications: 'close' ('C'), 
of which we find 45, and 'detached' ('D'), of which we find 37. For the purposes of this paper, 
'close' systems are defined as those for which there is no apparent distinction between in- 
and out-of-eclipse regions in the light curve, i.e., where there is no clear discontinuity in the 
slope of the curve. These should represent, for the most part, contact and over-contact binary 
systems. Those for which such a distinction is apparent are considered 'detached', though 
this category can also include semi-detached and near-contact systems. In some cases semi- 
detached or near-contact systems may also fall into the 'close' category instead, depending on 
the exact nature of the light curve. We have not attempted to further sub-classify the binary 
systems to avoid ambiguities that are better addressed with full light-curve modeling. 



Orbital parameters: Orbital periods, P, ar e obtained by calculating the Plavchan periodogram 
of the light curve (jPlavchan et al.ll2008l ). using a stand-alone version of the NASA/IPAC/ 
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NExScI Star and Exoplanet Database (NStED) periodogram servicel^ The light curve is 
folded on the period corresponding to the highest peak in the periodogram, and visually 
inspected to confirm the period is reasonable. If not, the next highest peaks are investigated 
until a plausible folding is found. The error in the period is estimated from the width of a 
Gaussian fit to the corresponding peak. In cases where the shape of the peak was such that 
a good Gaussian fit could not be made, the approximate peak was estimated by hand, and 
an appropriately large error assigned. 

The time of mid-eclipse, Tq, for the systems is measured by fitting a symmetrical inverted 
trapezium to the deepest and best-sampled minimum of each light curve and calculating the 
center time of the fitted trapezium. This fit also yields the formal errors for the eclipse center 
time. 

Distance: W UMa systems display a relatively well-established empirical period-color-luminosity 
relationship resulting from their approximately uniform-temperature common envelope: Ke- 
pler's law constrains the orbital radius as a function of period, and the orbital radius in turn 
constrains t he radiating su rface area and, hence, the luminosity and absolute m agnitude of 



the sy stem (jRucinskil l2004l ) . They therefore can be used as distance estimators. lEker et al 
(I2OO9I ) provided an updated calibration of the relationship based on 2MASS {J — H and 
H — Kg) colors; we use this relationship to estimate the distances for the 'close' binaries in 
the table where the detections have 2MASS counterparts, on the assumption that the vast 
majority are W UMa systems. 



Bilir et al.l (120081 ) provided a similar empirical relationship between 2MASS colors and abso- 
lute magnitude for main-sequence detached binaries, analogous to the color-luminosity rela- 
tionship for single main-sequence stars. We apply this relationship to the 'detached' binaries 
in the table, on the assumption of dwarf star status. 

Errors in the derived distances are estimated by standard propagation of the 2MASS photo- 
metric errors through the respective empirical relationships. The errors in the coefficients in 
the relationships are neglected, but the reported standard deviations about the relationships 
(0.26 mag and 0.49 mag for the W UMa and detached-system calibrations respectively) are 
added in quadrature to account for astrophysical effects that are not included in the relation- 
ships. In both cases, the effects of reddening are also n eglected: for our ass umed value of 



Ay at the distance of 25 Ori, using the extinction laws of lCardelli et al.l (|l989l ) gives redden- 
ing values of E{J — H) = 0.027 mag and E{H — K) = 0.022 mag, which in most cases are 
comparable to or less than the photometric measurement errors. Neglecting reddening does 
lead to a small systematic underestimation of the distances, generally by of order half the 
measurement error in most cases for our assumed mean Ay, and always less than the full size 
of the error, but it is simpler than attempting to realistically estimate the reddening for each 
individual source. 



^http: //nsted. ipac . caltech. edu/periodogram/cgi-bin/Periodogram/nph-simpleupload| 
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We note that these distance relationships are hkely to be inaccurate at PMS ages. Low 
mass stars (< 1.5 Mq) younger than lOMyr wih stih be on PMS tracks and can be expected 
to be brighter than predicted (albeit the close companions in our sample may complicate 
the evolutionary picture). The agreement between the list of young candidates from the 
color-magnitude diagram and the list from the distance estimates, however, provides some 
suggestion that the estimates are still meaningful in these cases (see sections 15.2.11 and I5.2.2P . 
Some systems also re present extrapolat ions beyond the forma l calibrated limits of the rela- 
tionships provided by lEker et al.l ((20091) and iBilir et al.l (|2008l ) , particularly for the detached 
systems at apparent types later than mid-K. The distances provided in the table therefore, 
are only intended as a guideline and should be treated with appropriate caution. 

R magnitudes: A single representative value for the magnitude, -Rmedj is given as a straight 
median of the measured differential R magnitude for each light curve. We estimate this to 
be accurate on an absolute scale to « 0.03-0.05 mag for chips 0-4 and 6-10, and ^ 0.07 mag 
for chips 5 and 11 (see section W72\i . The light-curve amplitude variation, AR, is calculated 
as the difference between the 1st and 99th percentiles of the measured R magnitudes for each 
light curve, in order to reject possible unflagged outliers and allow for a c ertain amount of 
statistical deviation. This is similar to the range used bv lBasri et al.l (|2010l ). 
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Fig. 4. — : Folded light curves for the eclipsing binaries identified in the data (section [5.1l and table [T]). The 
gray scale represents time of observation to give a sense of any long-term variation, running from 
black at the beginning of the PTF Orion observing run to pale gray at the end. Error bars are 
omitted for clarity. In cases where insufficient coverage was obtained to determine a period, or 
where a regular period could not be found, data covering the best eclipse detection are instead 
plotted in black directly against heliocentric Julian date (complete light curves for these are 
shown in figures I MTU)) . See tabled] for the full precision in the measured periods. (Continued 
in figures [5HH1) 
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due to an insufficient number of eclipse detections or because the data did not exhibit a regular 
period. Continued in figure [101 
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Fig. 10. — : Unfolded light curves - continued. 5-6432 (bottom) shows sparse coverage because most of 
the data is saturated; the light curve drops below the saturation threshold during the eclipse, 
however. 



Table 1. New Eclipsing Binaries from the PTF Orion Data 
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To 
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± 





.02 


11, 


,89 


±0, 


,02 


D 


0, 


,679631 


± 


0, 


,001025 


205, 


,83808 


± 


0, 


,00025 


349 ± 83 


14 


,58 


0, 


,06 


2-1066^ 


81.49801 


1.96661 


J05255953+0157593 


13, 


,64 


± 





,03 


13, 


,04 


± 


0, 
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5-12446° 83.34683 


2.97482 


,105332318+0258281 


11, 


,99 


± 


0, 


,03 


11, 


,47 


± 


0, 


,02 


11, 


,34 


± 


0, 


,03 


C 


0, 


,236633 


± 


0, 


,000017 


204, 


,65544 


± 


0, 


,00003 


317 ± 41 


13, 


,46 


0, 


,39 


7-5291« 


80.80171 


1.35431 


J05231240+0121151 


12, 


.13 


± 


0, 


,02 


11, 


,44 


± 


0, 


,02 


11, 


.31 


± 


0, 


,02 


C 


3, 


,623108 


± 


0, 


,023858 


192, 


,86419 


± 


0, 


,00265 


1442 ± 183 


13, 


,95 


0, 


,11 


10-10597 


82.53469 


1.81344 


J05300832+0148489 


11, 


.77 


± 


0, 


,03 


11. 


,27 


± 


0, 


,03 


11, 


.13 


± 


0, 


,02 


c 


0, 


,299531 


± 


0, 


,000050 


192, 


,72915 


± 


0, 


,00008 


348 ± 45 


12, 


,96 


0, 


,18 


Other Binaries of Interest 


































































0-3542 


80.00204 


2.25436 


J05200049+0215159 


14, 


.27 


± 


0, 


,03 


13, 


,74 


± 


0, 


,02 


13, 


.51 


± 


0, 


,04 


c 


0, 


,284425 


± 


0, 


,000077 


199, 


,77639 


± 


0, 


,00023 


934 ± 122 


15, 


,62 


0, 


,28 


0-9653° 


80.29957 


2.93869 


J05211191-I-0256189 


14 


.70 


± 


0, 


,05 


14. 


,18 


± 


0, 


,04 


13, 


.88 


± 


0, 


,06 


c 


0, 


,910619 


± 


0. 


,001313 


202, 


,84000 


± 


0, 


,10000 


2273 ± 339 


16, 


,23 


0, 


,29 


9-2980' 


81.86439 


0.98081 


J05272748+0058506 


14, 


.39 


± 


0. 


,04 


14. 


,12 


± 


0. 


,04 


14, 


.11 


± 


0. 


,07 


D 


1, 


,190701 


± 


0. 


,003662 


175, 


,79416 


± 


0, 


,00062 


2923 ± 900 


14, 


,87 


0, 


,07 


11-7037 


82.87870 


1.50007 


J05313088+0130004 


15, 


.01 


± 


0. 


,05 


14. 


,50 


± 


0. 


,07 


14, 


.25 


± 


0. 


,09 


c 


0, 


,306948 


± 


0. 


,000228 


199, 


,74559 


± 


0, 


,00116 


1395 ± 235 


16, 


,09 


0, 


,13 


11-8774^ 


83.32166 


1.68053 


J05331720+0140498 


13 


.23 


± 


0, 


,03 


12. 


,62 


± 


0, 


,03 


12 


.47 


± 


0, 


,03 


C 


0, 


,215651 


± 


0, 


,000007 


180, 


,94621 


± 


0, 


,00004 


450 ± 60 


14, 


,83 


0, 


,78 



Other Binaries 



0- 


■1580 


79, 


,94125 


2, 


,02664 


J05194592+0201361 


14, 


,35 




0, 


,04 


13, 


,91 


± 


0, 


,04 


13,78 ± 0.05 


D 


2, 


,061758 


± 


0, 


,001051 


198, 


,80653 


± 


0, 


,00012 


1334 ± 376 


15, 


,56 


0, 


,52 


0- 


■4268 


80, 


,46153 


2, 


,33903 


J05215079+0220205 


13, 


,03 




0, 


,03 


12, 


,67 


± 


0, 


,02 


12.66 ±0.03 


D 


5, 


,167944 


± 


0, 


,176421 


202, 


,79638 


± 


0, 


,00003 


1233 ± 307 


13, 


,90 


0, 


,70 


0- 


■4407 


80, 


,17316 


2, 


,35519 


J05204154+0221186 


14, 


.60 


± 


0, 


,03 


14, 


.07 


± 


0, 


,04 


13.97 ±0.07 


D 


0, 


,701540 


± 


0, 


.001087 


189, 


,75110 


± 


0, 


,00039 


1296 ± 389 


16, 


,07 


0, 


,15 


0- 


■5197 


80, 


,45244 


2, 


,44664 


J052 14859+0226480 


16 


.20 


± 


0, 


,13 


15, 


.42 


± 


0. 


,11 


15.60 ±0.25 


C 


0, 


,273236 


± 





.000064 


198, 


,80522 


± 


0, 


,00047 


2006 ± 627 


17, 


,24 


0. 


,80 


0- 


■5673 


80, 


,06399 


2, 


,50293 


J05201536+0230104 


14 


.91 


± 


0, 


,05 


14, 


.39 


± 


0, 


,04 


14.34 ± 0.08 


D 


0, 


,543151 


± 


0, 


.000185 


175, 


,83572 


± 


0, 


,00037 


1786 ± 602 


15, 


,93 


0. 


,19 


0- 


■80368 


80, 


,03662 


2, 


,76290 


0927-0078105 
























D 


0, 


,597359 


± 





.000604 


204, 


,70458 


± 


0, 


,00048 




18, 


,09 


0, 


,75 


0- 


■8177 


80, 


,28231 


2, 


,77769 


J05210774+0246397 


15, 


,98 


± 


0, 


,08 


15, 


,82 


± 


0, 


14 


> 15.21 


C 


0, 


,304414 


± 


0, 


,000214 


204, 


,74716 


± 


0, 


,00047 




17 


,11 


0, 


,25 


0- 


■9790 


80, 


,42457 


2, 


,95301 


J05214195+0257109 


15, 


,96 


± 


0, 


10 


15, 


,77 


± 


0, 


14 


15,49 ±0.23 


C 


0, 


,326311 


± 


0, 


,000147 


199, 


,72763 


± 


0, 


,00036 


3842 ± 1112 


17 


,12 


0, 


,43 


1- 


■1223 


80, 


,94578 


2, 


,00642 


J05234700+0200225 


15 


,26 


± 


0, 


,05 


14, 


,80 


± 


0, 


,07 


14,87 ±0.13 


D 


1, 


,252181 


± 


0, 


,000334 


192, 


,80870 


± 


0, 


,00012 


3488 ± 1528 


15, 


,91 


1, 


,59 


1- 


■2659 


80, 


,80171 


2, 


,18552 


J05231240+0211079 


15, 


.99 


± 


0, 


,11 


15, 


.90 


± 


0, 


,19 


15.12 ±0.16 


D 


0, 


,960161 


± 


0, 


.002785 


189, 


,89204 


± 


0, 


,00025 


1045 ± 585 


16, 


,67 


0, 


,59 


1- 


■3638^ 


81, 


,08221 


2, 


,31406 


J05241975+02 18507 


15, 


.68 


± 


0, 


,07 


15, 


.46 


± 


0, 


,12 


15.43 ± 0.20 


D 












180, 


,78801 


± 


0, 


,00028 


5650 ± 3628 


16, 


,63 


0, 


,43 



Table 1 — Continued 
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Note. — A description of the columns is given in section 15.11 Corresponding light curves are shown in figures Bl llOl Distance estimates are omitted in cases where 
2MASS magnitudes reported are limits only. Where possible the matched 2MASS identifier is given based on the J2000 source coordinates; in cases where no 2MASS 
match was found, the corresponding USNO-Bl.O match running-numbcr-based ID is given instead, with the format 'NNNN-NNNNNNNN' (with no preceding 'J'). The 
full 2MASS identifier has the format '2MASS JHHMMSSss+DDMMSSs,' where HH, MM, and SSss represent hours, minutes, and seconds of right ascension, respectively; 
and DD, MM, and SSs represent degrees, minutes, and seconds of declination, respectively. Lower case 's' represents digits to the right of the decimal point. Pull PTF 
survey catalog IDs can be constructed based on these coordinates, using the format 'PTFl JHHMMSS.ss+DDMMSS.s.' Since precision astrometry was not the primary 
goal of the PTF Orion project, it is preferable to use 2MASS/USNO-B1.0 coordinates for this purpose rather than our PTF Orion measured coordinates. 

^Only one eclipse obtained - period indeterminate. 

''Partial coverage of only one (presumed) eclipse obtained - period indeterminate, Tq not well constrained. 



'^Tq not well determined due to incomplete coverage of eclipse. Measurement error may be unreliable. 
'^Only one (very good) eclipse observed; period based on out-of-eclipse variation. 

"Unusual light curve, possible semi-detached system. Listed distance is assuming a contact binary, but using a detached- system distance estimate yields a distance of 
268 ± 64pc - see section [5A3l 

'Very short period W UMa system - see section [5.4.41 

^Used center of secondary eclipse for Tq, owing to poor coverage of primary. 
'^Three clear eclipses, but unable to find coherent period - possible triple system? 
'Nearby second source in USNO-B; chance of slight contamination. 

J Secondary eclipse not evident, or primary and secondary eclipses indistinguishable - possible factor of two ambiguity in P. 
''Only two eclipses obtained - period ambiguous. 

'Apparent pulsating binary - short-period oscillations seen, with ~ 1/2 hr period. See section [5.4.41 

™Somewhat distorted light-curve shape - may actually represent stellar pulsations rather than a binary system (D. Bradstreet, priv. comm., 2010). 
"See section [4.21 for a discussion of zero-point accuracy. 
"Identified as a candidate variable star bv lKraus et al. 
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5.2. Identifying Young Binaries 

We identify young binary candidates in our sample with a color-magnitude selection, and then 
follow several approaches to check for consistency, detailed below. We find nine candidates which 
are worthy of followup (see beginning of tabled]). 



5.2.1. Color- Magnitude Selection 

In figure [11] we plot the complete set of 82 binary systems in a color- magnitude diagram, using 
our median measured R values and the single-epoch 2MASS Kg magnitudes (excepting those with 
poor Kg measurements, i.e., those with 2MASS quality flags worse than 'C or error bars > 0.1 mag 
- see section 12. 3p . Isochrones are over-plotted which are b ased on the ISiess et al.l (12000) models . 
These are tuned to fit the Pleiades main sequence at 100 Myr ( Guieu et ahlboiC : Stauffer et al.|[2007 : 
Jeffries et al.ll2007l). and then cori verted to R versus R — Kg using polynomial fits to young-star 
colors ()Kenvon &: Hartmannlll995l ) . The blue tri a ngles mark the locations of previously identified 
(single) WTTS stars identified by iBriceno et al.l (|2007l ) which are also found in our data, and in 
the 2MASS ca talog. They li e along the 10 Myr isochrone, in good agreement with the 10 Myr age 
determined bv lBriceno et al.l (|2007l ) also using Siess et al. isochrones, and supporting our placement 
of the isochrones in the diagram. 

The uncertainties of magnitudes and colors are dominated by the amplitude of the physical 
photometric variability. As a result the error bars in figure [11] are assigned to represent the overall 
variability of the light curves (half the peak-to-peak amplitude). In the absence of multi-epoch 
2MASS measurements, we assume the variability in Kg to be the same as that in the PTF R 
measurements for the purposes of the error bars, and since the R and Kg are not contemporaneous, 
we assume the errors to be independent. We treat the WTTS sources in the same way to allow for 
their variability as well. 

Those systems lying below the marked main sequence are taken to be unassociated with the 
25 Ori/OBla association. Those falling above the main sequence (the PMS region for the assumed 
25 Ori distance and reddening) are highlighted with star symbols, and identified as candidate young 
systems. This criterion yields nine systems: 0-7220, 2-1066, 2-1144, 4-5331, 5-12446, 7-5291, 7-7604, 
10-10597, and 11-5402 (see table[l]). For comparison, we ran a Besancon Galactic model simulation 
for the PTF Orion field using the default model parameter^ ( Robin et al. 2003 ). This allows us 
to estimate how many field stars might be expected to lie in the PMS region. Scaling the fraction 
of sources found to lie above the 330 pc main sequence in the simulation to the size of our binary 
sample yields ^ 2-5 expected field sources, depending on the exact location of the main sequence 
cut. 



^http : //model . obs-besancon . f r/ 1 



- 33 - 
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Fig. 11. — : Color- magnitude diagram for the binary systems identified in this paper ( circles/star symbols) , 
and for other sources previously identified as single WTTS sources by iBricefio et al.l (j200lT ) 
which are found in the PTF Orion data (triangles). Star symbols highlight candidate young 
Orion systems, labeled with their PTF Orion IDs. The dotted lines mark rough completeness 
limits; the solid line marks the main sequence at 330 pc, with the assumed mean reddening 
values for the 25 Ori association; the dashed lines represent PMS isochrones at 1, 3, and 10 Myr 
(see section [5'. 2 .ip . Two further sources, 2-822 and 5-5767, are also labeled; these are discussed 
in section [5.4.21 
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There is an ambiguity with regard to binary systems on the color-magnitude diagram in 
that the mass ratios of the systems are unknown. For systems with a small mass ratio, one star 
dominates the luminosity, and hence the system's location on the diagram would fall as expected 
for a single star of the same type and color (diverging evolutionary effects due to the binarity 
notwithstanding). However, for mass ratios near unity, the two stars contribute similar fluxes, 
appearing up to ~ 0.75 mag brighter in total than a single star of the same color. There is a 
possibility, therefore, that some binaries could be shifted systematically upward in the diagram 
compared to single stars, and that location near the OBla/25-Ori group sequence could simply 
be a result of field binaries suffering from this effect. These nine systems should therefore only be 
considered a list of good candidate young systems from the sample, for which further spectroscopic 
follow-up is needed for proper confirmation. At the same time, we note that shifting downward by 
0.75 mag any of those candidates on the low edge of the cluster sequence in the diagram (assuming 
that R — K remains the same) still places them near the main sequence for the distance of the 
association, and does not rule out Orion membership. It is possible, for example, that binarity has 
affected their evolution compared to that of single stars, particularly given the close separations 
implied by the relatively short orbital periods of our sample. In the following sections, we investigate 
other consistency tests for the proposed young ages of the candidates. 



5.2.2. Distances 

Figure [12] shows a histogram of the log of the derived distances for the binaries in table [J 
(excluding those with 2MASS quality flags worse than 'C ~ see section 12. 3p . The number of 
binaries per logarithmic bin increases with distance, peaking at around log(distance/pc) = 3.1 and 
then tailing off. This we attribute to the background field distribution of binaries in the galaxy in 
combination with the sensitivity limit of the sample and the logarithmic bin sizes. In addition we 
se e a distinct pe a k at around 330 pc, the distance to the OBla and 25 Ori associations adopted 



by lBriceho et al.l (|2005l . 120071 ). We take systems with derived distances falling at the 25 Ori/OBla 
distance within a factor of 1.5 times their respective errors, and with errors below 150 pc, to be 
good candidates for membership in the association. This criterion yields exactly the same list of 
candidates as the color-magnitude diagram selection, consistent with their proposed young nature, 
with the exception only of source 7-5291. This source presents some peculiarities, however, and 
there are indications that the initial distance estimate is probably incorrect. An estimate based 
on the detached-system calibration is likely more appropriate than the presumed W UMa distance 
calibration, and satisfies the above criteria. We therefore do not rule it out as a young candidate 
(see section r5. 4. 31 for more discussion). 

The same mass-ratio ambiguity as discussed in the previous section can also affect the distance 
estimates, although probably to a lesser extent since the distance calibrations used are based on 
broad samples of binary systems, and therefore should represent the relation for a 'typical' mass 
ratio. iBilir et al.l (j2008l ) report a standard deviation of 0.49 mag about the relation for detached 
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Fig. 12. — : Histogram of estimated distances to tlie binary systems (section [5221) ■ ^ peak is seen around 
330 pc, the estimated distance of the 25 Ori group and Orion OBla (logj^Q distance = 2.52, 
marked with a dashed hne). Bin size is 0.15 log^^g (distance/pc). 



systems, pointing out that this could be largely e xplained by the mass-ratio ambiguity, which 
should lead on average to an error less than 0.4 mag; Eker et al. (j2009l ) report a standard deviation 
of 0.26 mag for W UMa systems. In the worst possible case where a mass ratio of unity leads to 
a f» 0.75 mag underestimate of the distance modulus, this would cause a distance underestimate 
of ~ 40%. Even at this extreme, however, an estimated distanc e at the 330 pc ma, r k still places 



the t rue value near the more distant OBlb association (~ 440 pc: iBricefio et al. 
19991 ). and such a source would remain of potential interest. 



2005 



Brown et al 



5.2.3. Stellar Variability 

Intrinsic stellar variability can also be a usefu l indicator for you ng stars, which tend to be 
more spotted and active (e.g., IBricefio et al.l l200ll : iLamm et al.l l2004l ) . In figure [13] we plot an 
estimate of the intrinsic variability of the detached binary systems in our sample for which a 
meaningful estimate could be calculated. All 'close' systems are omitted, as is 5-961, a likely 
semi-detached system, since they show continuous strong ellipsoidal variation which confuses the 
intrinsic variability measurement. To quantify the intrinsic stellar variability, the out-of-eclipse rms 
R magnitude is first calculated as the outlier-resistant standard deviation for a given light curve 
between orbital phases 0.1-0.4 and 0.6-0.9, effectively cutting out the eclipses in the detached 
binaries (assuming circular orbits, which is reasonable given the short orbital periods in our sample). 
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To account for the apparent noise floor in figure[3]a constant 0.004 mag is added in quadrature to the 
measurement errors. Variabihty is then calculated by subtracting the median measurement error for 
each light curve (in quadrature) in order to leave only intrinsic variability. This quantity provides 
an estimate of the stellar activity, though with the caveat that low-level ellipsoidal variations may 
also be included in the metric. The value of the reduced compared to a constant R magnitude 
at the weighted mean for each source is also plotted, to provide an indication of the statistical 
significance of the measure. 

The effects of photon error having been removed from the variability diagram, the general 
increasing trend with fainter magnitudes is attributed to the tendency toward increased activity at 
later spectral types, since more of these are seen at the faint end of the sample. The same candidate 
young binary systems (section l5.2.ip are again marked with star symbols (excluding the three 'close' 
systems, 7-5291, 5-12446, and 10-10597). They are intrinsically more variable on average than the 
other detached systems in the sample, consistent with their suspected young ages. (The elevated 
source below 11-5402 is an additional system that shows marked variation probably caused by hot 
or cold spots (7-337 - see figure [6]).) 



5.2.4- Proper Motions 



Proper motions provide a further consistency check on the candidate young binary systems, 
although for the 25 Ori cluster and the broader Orion region, the proper motion is near zero, making 
it hard to distinguish b etween 25 Ori association members, Ori QBl members, and more distant 
background sources (see lde Zeeuw et al.lll999l : lKharchenko et al.ll2005l ). Figure [14] shows the proper 
motions of each of th e binary systems, ta ke n from the recently releas ed PPMXL catalog of positions 
and proper motions (jRoeser et al.ll20ld ). iKharchenko et al.l (120051 ) reported a proper motion for 
the 25 Ori cluster (identified as ASCC 16 in their catalog) of 0.75 it 0.22 mas yr~^ in r ight ascension 
and —0.18 zt 0.29 mas yr~^ in declination based on the brighter stars in the cluster. |Briceno_et_al 



( 20071 ) report a ibl.7kms~^ width to the core of the peak of their distribution of radial velocity 
(RV) measurements for the cluster, which corresponds to a ±1.09 mas yr~^ dispersion in proper 
motions at 330 pc (indicated in the figure), assuming a three-dimensionally isotropic distribution 
of velocities. 

The nine young candidate systems are again highlighted with red star symbols. Their proper 
motions are largely consistent with 25 Ori/OBla membership, with the exception of 5-12446, with a 
proper motion of 30.0 it 3.8 mas yr~^ in right ascension despite its estimated distance of 317 it 41 pc, 
which is extremely close to that of the 25 Ori/OBla association. This suggests it is probably an 
interloper rather than a true member (see section r5.4.3p . 
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Fig. 13. — : Intrinsic stellar variability estimates for all the binary systems classified as detached (section 
15.2. 3p . Star symbols highlight those which are the candidate young systems selected from the 
color-magnitude diagram (figure [11]) . Variability (top panel) is calculated as the quadrature 
difference between the out-of-eclipse rms and the median measurement-error for each light curve 
(with an additional estimated 0.004 mag noise floor added in quadrature); sources where the 
rms is below the median measurement error are omitted. The reduced against the weighted 
mean R magnitude is shown for the same data (bottom panel), including those sources with 
rms smaller than the measurement error. System 5-961 is omitted from both panels since it 
shows significant ellipsoidal variation (figure [5]) . 
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Fig. 14. — : PPMXL proper motions for all the identified eclipsing binary star systems. Candidate young 
systems are marked with star symbols, and have error bars indicated. Error bars for all other 
systems are omitted for clarity, but are similar in size. The squ are symbol marks the proper 
motion of the 25 Ori cluster measured bv lKharchenko et al. ( 20051). with error bars representing 
the dispersion derived from RV measurements by Briceno et al.l ( 2007 ). (See section 15.2.41 ) 



5.2.5. Infrared Analysis 

Finally, we also searched for any existing Spitzer IRAC/MIPS data, to check for any infrared 
excess that might be indicative of T-Tauri-like disks in the nine candidate young systems. Data 
were found for 2-1066, 2-1144, 7-5291, 7-7604, and 10-10597. These are listed in table El and SEDs 
are shown in figure [T5| where the gray lines represent an extrapolation of the Rayleigh- Jeans slope 
from the Kg measurements to guide the eye with regard to the expected photospheric flux. A 
color difference between any two points on such a gradient will have a value of zero, and anything 
lying above the slope is suggestive of an IR excess. As discussed in section 12.31 comparing to 
Ks is not optimal, but [3.6] is in many cases not available. 2-1066 seems to show a clear excess, 
with Kg — [24] = 4.75 ± 0.06 mag. Though it lacks a [3.6] measurement, the two IRAC bands 
which are measured ([4.5] and [8]) also lie above the extrapolated photosphere as can be seen in 
the figure, supporting this conclusion (though it is also possible that the excess may be caused by 
confusion from the secondary in the system). The remainder show no detected excess, or at most, 
the possibility of a marginal very weak excess. These results are not surprising: by design, the age 
of the region in the field is such that for the most part the disks have largely dissipated. 



5.3. Identifying Low-mass Binaries 

In figure [TBI we show a JHKg color-color diagram plotting all of the 82 eclipsing binary systems 
for which good quality photometry are found in the 2MASS point source catalog. Also indicated 
are the regions encompassed by the main sequence and the giant branch. The nine young binary 



Table 2. Existing Spitzer Data for Candidate Young Binaries 



ID 


J 


H 


Ks 


[3.6] 


[4.5] 


[5.8] 


[8] 


[24] 


Ks - [24] 


[3.6] - [24] 


2-1066 


13.64±0.03 


13.04±0.03 


12.76±0.03 




11.84±0.06 




11.24±0.06 


8.01 ±0.05 


4.75±0.06 




2-1144 


13.09±0.03 


12.43±0.03 


12.20±0.03 




11.82±0.06 




11.77±0.06 


11.68±0.34 


0.52±0.34 




7-5291 


12.13±0.02 


11.44±0.02 


11.31±0.02 


11.20±0.06 




11.21±0.06 




10.47±0.12 


0.84±0.12 


0.73±0.13 


7-7604 


13.09±0.03 


12.39±0.02 


12.18±0.03 


11.97±0.06 


11.87±0.06 


11.87±0.06 


11.89±0.06 


10.18i0.ll 


2.00±0.11 


1.79±0.13 


10-10597 


11.77±0.03 


11.27±0.03 


11.14±0.02 


10.93±0.06 


10.94±0.06 


10.92±0.06 


10.92±0.06 


11.48±0.20 


-0.34±0.20 


-0.55±0.21 



Note. — All candidate young binaries are listed for which Spiteer IRAC/MIPS data could be found (scction l5.2.5ll . 2MASS J,H, and Ks measurements 
are repeated from table [l] for comparison. All values are in magnitudes. 
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Fig. 15. — : SEDs for the candidate young binary systems for which IR Spitzer data could be found (section 
15.2.51 and table [2|). The SEDs are all in log XF\ in cgs units (erg s^^ cm^^), against log A in 
microns. Crosses are median PTE R, diamonds are 2MASS {JHKs, 1.25, 1.65, and 2.17 /im), 
circles are IRAC (3.6, 4.5, 5.8, 8/im), and squares are MIPS (24 ^m). A gray line is extended 
from Kg to guide the eye, representing an extrapolated Rayleigh- Jeans slope. If the object is 
a star, there is no excess at Ks, no significant reddening, and the star is earlier than mid-M, 
this would be the approximate location of the photosphere. 2-1066 shows a clear possible IR 
excess; a very weak excess appears to be present at 24 /xm for 7-7604. 



candidates are highlighted with star symbols as before. 

From the figure we can see that, with the exception of the two contact binaries, 5-12446 and 
10-10597, all of the previously discussed candidate young systems also have 2MASS colors that - 
assuming no major reddening - are consistent with M dwarf primaries. These systems are therefore 
consistent with being both young and low mass. 

In addition to the candidate young systems already discussed, we also mark on the diagram 
the five other systems which fit in or near the M-dwarf region of the plot, but are not necessarily 
young or associated with Orion: 0-3542, 0-9653, 2-822, 5-5767, and 11-7037, of which three (0-3542, 
0-9653, and 11-7037) are classified as 'close' systems. These five systems are also broken out at the 
beginning of table [H and are discussed in sections 15.4.21 and I5.4.4I 

Dereddening any of these systems by the mean assumed reddening, indicated on the plot, 
would make little difference to their status. For the young-system candidates, it only brings the 
sources closer to the main sequence. However, it remains possible that some sources are either 
obscured by tight localized regions of interstellar dust, or intrinsically reddened objects. The small 
IR excesses in the Spitzer data for 2-1066 (section [5.2.5p . and possibly 7-7604, may be indicative of 
intrinsic reddening. It should be noted that the possibility of unresolved source confusion remains, 
which could also lead to a distortion of the measured colors. 
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Fig. 16. — : Color-color diagram for the new binaries based on 2MASS colors (section I5.3[) . The green 
hatched region marks the main sequence; the blue region marks the giant branch; and the red 
marks the beginning of the L dwarf regime. The assumed mean reddening vector for the field is 
indicated. The candidate young binary systems are marked with a star symbol as in previous 
figures. In addition, five other potential low mass binary systems are marked with triangles on 
the basis of their location in the diagram. 
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5.4. Discussion 



In the previous sections ()5.2I and I5.3P we discussed the approaches employed for identifying 
particular sources of interest in the eclipsing binary sample. We identified nine PMS candidate 
systems (0-7220, 2-1066, 2-1144, 4-5331, 5-12446, 7-5291, 7-7604, 10-10597, and 11-5402), of which 
two (5-12446 and 10-10597) appear to be contact binary systems, one is also categorized as 'close' 
but is somewhat anomalous (7-5291), and the remainder ar e likely low-mass systems . (Co ntact 



binaries are not normally found as late as M-type - see, e.g.. lRucinskill2004l : iBilir et al.ll2005l ). We 
also identified three notable very red contact-binaries (0-3542, 0-9653, and 11-7037). In addition, 
two more systems (not previously discussed) are particularly interesting: 11-8774, a contact binary 
with an exceptionally short period, and 9-2980, a detached system which appears to exhibit stellar 
pulsations. We summarize the results here and discuss some of the specific systems. Since there is 
significant overlap between the results from the analysis above, we break down the specific sources 
identified into the four categories below, which are reflected at the beginning of table [TJ 



5. 4-1- Young Low-mass Binary Candidates 

Six of the nine young candidate systems we found fell also into the category of candidate low 
mass (MO or later) systems: 0-7220, 2-1066, 2-1144, 4-5331, 7-7604, and 11-5402. These specific 
systems are discussed in more detail below. The remaining three are close binaries, discussed in 
section [5.4.31 (though 7-5291 may also qualify as low mass). 

0-7220: Presents an unusual light curve. The regular variation outside the primary and (small, 
~ 1%) secondary eclipses is comparable with the depth of the primary eclipse itself (see figure 
2]), and its phase is not consistent with ellipsoidal variation. It may be attributable to star 
spots: either one spot system with the star rotating at twice the orbital period, or more likely, 
two spot systems, separated by ^ 180° in longitude, rotating synchronously with the 0.680 
day orbit. This would also be consistent with the long-term variation of the out-of-eclipse 
modulation during the length of the observations, which can be seen from the gray scale in 
the light curve plot (figured]). 

2-1066: Although coverage for this object was good, only one eclipse was observed (figuredl bottom 
right), so the period is indeterminate. The eclipse is flat-bottomed, making it a noteworthy 
source for further observations. Out-of-eclipse variation is on the level of ~ 4%. It is also 
notable for its signiflcant apparent IR excess, with Kg — [24] = 4.75 ± 0.06 mag (section 15.2.51 
and figure [T5]) , which is suggestive of a disk (although it may also be attributable to the 
secondary companion). 

2-1144: Displays a similar light curve to 0-7220, though a secondary eclipse is not clearly evident 
in the data, and the out-of-eclipse variation peaks only once per orbit. Again, this is not 
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consistent with ellipsoidal variation, and is likely attributable to star spot/hot spot variation 
with the stellar rotation period synchronous with the 0.554 day orbital period. 

4-5331: Only part of one eclipse was captured for this object (figure [5|), and therefore the period 
is not constrained, although coverage was good. The rest of the light curve is variable with 
an amplitude of ~ 0.15 mag, with no clear period (figure [9]). It is particularly noticeable for 
its very red colors, located as late as M9 in figure [161 it is likely significantly reddened since 
if it had such a low mass it would appear much fainter at the distance of Orion. No coverage 
was found in the Spitzer data, but this reddening could be intrinsic, perhaps due to a disk. 
Though it seems likely that a stellar eclipse is the best explanation, it remains possible that 
the sudden drop in brightness is not in fact an eclipse from a binary companion, but that 
4-5331 is a TTS, perhaps a 'dipper'-like CTTS similar to 9-6886, shown in figure [Ml (section 

[SD- 

7-7604: Shows a small but significant level of flaring consistent with its proposed young age and 
a low mass - hence the apparent discontinuity in the folded light curve (figure [BJ. The 
object shows a possible very weak IR excess in the Spitzer data at Kg only, with Kg — [24] = 
2.00±0.11mag and [3.6] - [24] = 1.79±0.13mag (table[2land figure [T5l). 7-7604 lies closer to 
the giant branch than the main sequence in the JHKg color-color diagram (figure [T6l) . but if 
it is intrinsically reddened, this would bring its true color closer to the ~M0 main sequence 
region. The apparent broadness of the eclipses suggests quite a closely separated system, with 
significant tidal distortion. Similarly to 7-5291 (see below, section [5.4.3p . its short 2.07 day 
period argues against it being a giant system, since a giant would be larger than the implied 
orbital separation. 

11-5402: Shows clear, roughly sinusoidal out-of-eclipse variation (figure [HI), again most likely due 
to star spots given the phase of the variation relative to the eclipse. Only one eclipse was 
observed for the object. The eclipse has good coverage despite being narrow, and is shown 
in more detail in figure ITTl but the period of the binary is therefore not unambiguously 
determined. We take the period of the out-of-eclipse variations as the best guess of the 
orbital period, on the assumption of star spots co-rotating with the orbital period as for 
0-7220 and 2-1144 (though the period here is longer, at 5.5 days). 



5.4-2. Other Low-mass Binary Candidates 

We identified two other detached binary systems as low-mass candidates on the basis of the 
2MASS JHKg color-color selection (figure [T6l section [531) • Given their distance estimates, and the 
fact that they do not appear in the PMS region of the color- magnitude diagram (figure ITTl) , they 
are most likely background systems unconnected with the Orion association. 

2-822: Shows a straightforward grazing-transit detached light curve, with out-of-eclipse variation 
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Fig. 17. — : Stretched view of the Hght curve for 11-5402, showing the single echpse observed for the object. 

The period for this system is estimated from the periodic out-of-ecHpse variation (see section 
I5.4.1|) . Error bars are omitted for clarity. 

suggesting longitudinal temperature variation on one of the components, synchronous with 
the 0.9204 days orbital period (figure Hj). The JHKg colors suggest a spectral type of ssMSV 
(though see the caveat below). Dereddening by the assumed mean reddening vector makes 
little difference to this estimation; dereddening by a greater amount would move the system 
off the main sequence, unless it is a late F-dwarf system, which would require a rather extreme 
reddening {Ay ~ 5 mag). 

5-5767: Only two grazing eclipses were detected for this detached system, leaving a 'factors-of- 
two' ambiguity in the period, which is resolved by the periodic out-of-eclipse variation, if we 
assume that it is synchronous with the orbital period. Folding on the resulting 2.105 day 
period suggests that both are primary (or both secondary) eclipses. From the phase coverage 
(figure [5]) we see that the other eclipse was missed, assuming it would have fallen at phase 0.5. 
The out-of-eclipse variation appears to be caused by star spots. The JHKs colors suggest a 
spectral type of RiM3V, unless it is a heavily reddened G-type system, with Ay ~ 2.5 mag. 
Again, see the caveat below, however. 

In both cases, we note that given the sizes of the their respective 2MASS measurement errors, 
the spectral type estimates are quite uncertain. The same two systems are labeled in the color- 
magnitude diagram in figure [TTJ At face value, their R — Kg values suggest somewhat earlier 
spectral types, around late K to MO depending on reddening. However, the difference in epochs 
for the R and Kg observations leads also to large errors in ii — K^, as reflected in the error bars; 
given the errors in both the color-magnitude and the color-color diagrams, the measurements are 
not inconsistent. The true spectral types likely lie somewhere between the two values, around early 
M type. Followup spectroscopy is needed to provide a better estimate. 
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5.4-3. Young 'Close' Binary Candidates 

The majority of the 'close'-type binary systems in our sample appear most likely to be contact 
binaries from the shape of their light curves. Two W UMa systems are identified in our binary 
sample as candidate young 25 Ori/OBla-association members. It remains possible that their colors 
are distorted by unresolved third sources within the PSF, which could confuse their apparent 
young status, but if their ages are in indeed in the 7 — 10 Myr range, then there would be significant 
constraints on the mechanisms by which they could have formed, and they may merit further 
investigation. A third candidate young 'close' system is less likely to be a contact binary, but 
shows some unusual properties. 



5-12446: Shows a straightforward W UMa-type light curve (figure [6]), with period 0.23663 days. 
This system is an outlier in proper-motion space, with a proper motion of 30.0 it 3.8 mas yr^^ 
in right ascension (figure [T^ . It is probably an interloper; given its location somewhat below 
the cluster sequence in figure [TT\ it could also be an unrelated system with a mass ratio 
near unity that has been shifted upward in the color magnitude diagram due to the resulting 
luminosity enhancement, as discussed in section [5. 2. 1[ The expected location of such contact- 
systems at young ages with respect to single-star isochrones is uncertain, however, and the 
system's estimated distance (317 ± 41 pc) is extremely close to that of the 25 Ori cluster. 
Even if this is an underestimate due to a near- unity mass ratio (section I5.2.2p , it may still 
lie in the range of the bro ader Orion a s sociat ion. The system was previously identified as a 



candidate variable star bv lKraus et al.l (j2007l ). 



7-5291: This object presents a peculiar light curve. Classified as a 'close' type system under 
our definition, it has the appearance of a contact binary system, but with an unusual form. 
Assuming the two minima are indeed eclipses, then the width of the eclipses alone would 
suggest anything from a contact to a semi-detached system. However, the maxima and 
minima are rather sharp, and the primary ingress and egress are somewhat asymmetric. The 
period of the system, 3.62 day, is much longer than that expected for a W UMa system, and 
appears as an outlier in figure [181 The 2MASS JHKg colors are consistent with a mid- to 
late-K-type giant or possibly a late K or MO dwarf primary (see figure [TBI) . but the period 
does not allow for either of these: a dwarf system would not even be significantly tidally 
distorted at such an orbital separation (~ lOi?0), and therefore could not be near contact; a 
giant would be larger than the orbital separation. This suggests a subgiant-sized primary. 

One possibility is that distorted shape of the light curve results from a semi-detached (ac- 
creting) system with a disk. A semi-detached system would also be more consistent with 
the marked difference between the two eclipse depths, indicative of a relatively large tem- 
perature difference. If we assume that the system cannot be a contact binary owing to its 
period, then using the alternative (detached) distance calibration puts it at a distance of 
268 ± 64 pc, which is ~ Icr from the 25 Ori/OBla distance. A lOMyr-old infiated young 
low-mass primary, however, still would not be large enough. An alternative explanation for 
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the light curve may be periodic occultation by a warp ed disk, in an analogue to AA Tau 



(jBouvier et al.ll2003l ). similar to the objects reported by iMorales-Calderon et al.l (j201ll ). The 



available (single-epoch) S'piteer photometry of this object (figure [151 section [5.2.5p . however, 
shows at most only a very marginal IR excess at 24 //m, with Ks — [24] = 0.84 it 0.12 and 
[3.6] — [24] = 0.73 ±0.13 mag (table[2]). The estimated spectral type is on the border between 
K and M, however, and source confusion at the relatively low spatial resolution of the MIPS 
data is a potential concern. Additional followup is needed to better constrain the spectral 
type and Ay and better understand this object. 

10-10597: Again shows a straightforward W UMa-type light curve (figure [8]), with a period of 
0.29953 days. Spitzer data were found for this object showing no IR excess; this is unsurpris- 
ing, as we might expect any disk to be disrupted in such a system. Like 5-12446, this system 
also lies somewhat below the cluster sequence in figure [TTl with a similar distance estimate 
(348 lb 45 pc), so the same comments with regard to possible near-unity mass ratio may apply. 



5. 4-4- Other Binaries of Interest 

Among the remaining binaries, we identified a further five which are not necessarily young or 
associated with Orion, but which are noteworthy for other reasons. 



0-3542, 0-9653, and 11-7037: In the JHKs color-color diagram (figure [H]), while all the other 
low-mass (MOV and later) systems are classified as 'detached' and lie on or very close to the 
main sequence, these three alone lie somewhat below the M-dwarf branch, and show W UMa- 
like light curves. If they are genuine contact binaries and not near-contact or semi-detached 
systems, their location on the diagram suggests that they are either anomal ously low-mass 



systems (W UM a systems are not normally found as late as M-type - see, e.g. lRucinskill2004l : 



Bilir et al.ll2005l ). or probably more likely, particularly reddened W UMa systems {Ay ~ 2- 



3 mag) of earlier spectral type. In either case they may be interesting for further investigation. 
Given thei r unus ual colors, we note that they lie outside the formal validity range of the 



Eker et al.l (|2009l ) calibrations, and their respective distance estimates may be u nreliable 



0-9653 was previously identified as a candidate variable star bv lKraus et al.l (|2007l ) 



9-2980: Shows small but pronounced ~ 0.01 mag, ~ 45min modulations superimposed on the 
eclipsing binary light curve. Stellar pulsations seem the most likely cause. The 2MASS colors 
of the source imply a spectral type of ~G0, and the significant out-of-eclipse variations suggest 
both spot modulation and ellipsoidal variation indicative of a tidally distorted system. More 
observations are needed to sample the secondary eclipse. 

11-8774: Figure [T8l shows a histogram of the periods of all the binaries from table [H highlighting 
the distribution of the 'close' type systems. W UMa systems have a known sharp cutoff in the 
period distribution at the short-period end below around 0.22 days, just below a peak in the 
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Fig. 18. — : Orbital period Iristogram for all eclipsing binaries for which a period estimate could be found 
(section 15.4. 4|) . The solid regions indicate the fraction which are 'close' (probably mostly 
contact binary) sources; the hatched regions indicate the remaining detached sources. At the 
shortest end of the distribution is 11-8774, with among the shortest W UMa periods known 
(section [5.4.41) . The outlying 'close' source with a long period at just above log(period) — 0.5 
is due to object 7-5291, which presents a somewhat unusual light curve (see section [5.4. 3p . 



distribution (see, e.g.. iRucinskil 120071 ) . The cause of this cutoff remains largely unclear. We 
have not attempted to correct for the effects of the observing windowing function or magnitude 
completeness effects, but it is still clear that a significant number of our targets lie near the 
short-period end of this distribution, which is not currently well sampled (|Rucinskill2007l ) . In 
particular, the shortest-period source, 11-8774, has a period of 0.2156509 ± 0.0000071 days, 
putting it among the most rapidly orbiting W UMa systems known. 11-8774 also displays 
clear flat-bottomed primary eclipses, meaning the eclipses are total rather than grazing (figure 
[T9l) . This makes the object a particularly good target for constraining parameters via light 
curve modeling. The derived distance of 11-8774, 450 it 60 pc, may be too distant for the 25 
Ori cluster and the larger O rion OBla associati o n (?^ 330 pc), but is consistent with that of 
the more distant (~ 440 pc; iBricefio et al.ll2005l : iBrown et al.lll999l ) - and younger - OBlb 
association. It could perhaps be an outlier of the OBlb association, although at such a 
young age, it would be far from the main-sequence, and the distance estimate would likely 
be systematically wrong. If, however, the system does lie within the general distance of the 
broader Orion association, it could constitute a young W UMa system. As such it would be 
useful in understanding and constraining the formation and evolution of such systems. 
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Fig. 19. — : Stretched view of ttie folded liglit curve for tlie very sliort period W UMa system, 11-8774, 
centered on the primary eclipse (see section [5.4.41 ). The flat bottom of the eclipse can clearly 
be seen. The gray-scale is again used to represents time of observation, running from dark 
gray at the beginning of the PTF Orion run to pale gray at the end of the run. The secondary 
eclipse (not shown) is also flat-bottomed. Error bars are omitted for clarity. 
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6. CTTS Sources 



Another class of light curve found in our data is that of CTTSs. These young systems exhibit 
particularly large and irregular brightness modulations, thought to be related to photospheric spot s, 
short-lived accretion hot-spots, and/or obscuration due to circumstellar dust (jHerbst et al.lll994l ). 
On the basis of this variability, we identified visually 16 candidate CTTS light-curves which showed 
irregular variability on scales of 1 mag or more, listed in table [31 These include 14 candidates 
which are n ew to our knowledge, alon g with one known CTTS, CVSO 35 (8-10691), previously 



reported by iBricefio et alJ (12005 



20071') . and one previously reported as a candidate WTTS, SDSS 



J052700.12+010136.8 (9-3485) (IMcGehed \20m . The light curves for all 16 of the sources are 
shown in figures [20]| 251 

Figure [26] shows the location of the candidate CTTS sources on the same color-magnitude 
diagram as previously shown for the binaries (figure [TT]) . All clearly lie in the PMS region for the 
distance and mean reddening of the 25 Ori region. 

Of the 16 CTTS sources, 10 were found with coverage in the Spitzer archive data (including 
the one known CTTS), and all show IR excesses suggestive of disk presence. The SEDs of these 
sources are shown in figure [271 All the objects appear to have clear excesses, and for those where 
multi-band Spitzer data are available, the excess is seen in more than one band, lending credence 
to the identification of a disk excess. For the three where there are only 24 fim detections, the 
apparent 24 /xm excess is very large {Kg — [24] = 3.52 it 0.06, 5.45 ± 0.06, and 4.69 ± 0.06, for sources 
5-9166, 5-9411, and 7-3471, respectively), and seems to be real, although the uncertainty is larger 
than the formal error on the points because of the uncertainty in Ay and spectral type. 

1-3452: Seems to show a slight discontinuity between JHKg and the Spitzer IRAC bands 
(figure [27]) . This may either be a result of the contribution from the black-body function of the 
disk, or due to intrinsic variability between the epochs of 2MASS and Spitzer observations. 

9-6886: Displays a light curve which appears flat for a substantial portion of the time, but with 
downward d eviations of up to almost a mag nitude (figure[23]). This is reminisc ent of 'dipper' object s 
discussed bv lMorales-Calderon et al.[ (|201l[ ). for which AA Tau is a prototype (jBouvier et al.|[2003l ). 
(These objects were previously alluded to in section [5.4.31 with regard to object 7-5291.) 

The new sources are spatially distributed fairly evenly across the PTF field of view, with no 
obvious clustering. This suggests that they are likely associated with the broader Orion OB la 
association rather than with the 25 Ori group itself, which is at approximately the same distance. 



Table 3. Candidate CTTS: Sources with Irregular Light Curves and AR ~ 1.0 mag or More. 



ID 



a S 
(deg) (deg) 



2MASS 



PM a PUS Kmed'* AR an J H Ks 

(mas yr~^) (mas yr~^) (mag) (mag) (mag) (mag) (mag) (mag) 

[3.6] [4.5] [5.8] [8] [24] Ks - [24] [3.6] - [24] 

(mag) (niag) (mag) (mag) (mag) (mag) (mag) 



1-3452 81.07590 2.28890 05241824+0217200 7.0±3.9 2.7±3.9 17.33 1.08 0.23 13.60±0.02 12.59±0.02 11.88±0.03 

10.36±0.05 9.86±0.05 9.56±0.05 9.17±0.05 6.91±0.05 4.96±0.05 3.45±0.07 



2-8756 81.42163 2.81835 05254119+0249059 0.3+3.9 -7.7+3.9 16.33 



1.09 



0.27 13.59+0.03 12.86+0.03 12.35+0.03 



4-2045 82.63182 2.05147 05303164+0203051 3.8+4.0 -0.3+4.0 

9.86+0.05 9.58+0.05 



14.49 1.13 0.30 11.97+0.02 

9.39+0.05 8.62+0.05 5.29+0.05 5.19+0.05 



11.02+0.02 10.48+0.02 
4.57+0.07 



5-9166 82.87327 2.68835 05312959+0241183 -0.8+4.0 0.2+4.0 15.88 



0.92 0.24 13.19+0.03 

• • • 8.35+0.05 3.52+0.06 



12.39+0.03 11.87+0.03 



5-9411 82.97474 2.70859 05315396+0242310 6.1+3.8 -4.4+3.8 17.55 



1.72 0.45 14.97+0.05 

7.73+0.05 5.45+0.06 



13.99+0.04 13.18+0.04 



7-3471 80.97543 1.13702 05235410+0108128 0.5+3.9 -10.9+3.9 16.03 



1.14 0.20 13.50+0.03 

7.60+0.05 4.69+0.06 



12.76+0.03 12.29+0.03 



8-6391 81.26950 1.34760 05250468+0120509 -0.4+3.9 

10.06+0.05 



-1.6+3.9 



15.37 
9.21+0.05 



0.90 0.24 12.46+0.03 

5.15+0.04 5.91+0.05 



11.67+0.02 11.06+0.02 
4.92+0.07 



8-10691" 81.44118 1.76381 05254589+0145500 6.1+3.9 -1.7+3.9 

9.74+0.05 9.42+0.05 



13.77 1.60 0.53 11.46+0.03 

8.94+0.05 7.58+0.05 4.73+0.04 5.62+0.05 



10.76+0.03 10.35+0.02 
5.01+0.07 



9-3150 82.23882 0.99533 05285732+0059431 5.6+16.3 -31.4+16.3 18.70 1.38 0.30 13.90+0.03 13.26+0.03 12.89+0.03 



9-3485'' 81.75044 1.02696 05270012+0101368 2.3+3.9 -3.2+3.9 16.46 



0.78 



0.18 13.18+0.02 12.47+0.02 12.10+0.03 



9-6886"= 82.08428 1.35453 05282024+0121159 3.5+3.9 -0.4+3.9 15.07 0.71 0.22 13.05+0.02 12.31+0.02 12.06+0.02 



Table 3 — Continued 



ID « <5 2MASS PM a PM <5 -Rmod'^ AR an J H Ks 

(deg) (deg) (mas yr^-^) (mas yr^-"-) (mag) (mag) (mag) (mag) (mag) (mag) 

W^] W\ ^] Ks - [24] [3.6] - [24] 

(mag) (mag) (mag) (mag) (mag) (mag) (mag) 



9- 10304 81.73075 1.67283 05265537+0140224 -16.0±3.9 2.9±3.9 15.16 0.82 0.13 12.85±0.03 11.96±0.03 11.42±0.02 

10.59±0.05 10.34±0.06 10.24±0.06 9.76±0.06 6.08±0.05 5.34±0.05 4.51±0.07 

10- 2160 82.72717 0.94125 05305452-1-0056286 5.4±4.0 -1.5±4.0 17.10 0.76 0.14 13.34±0.03 12.80±0.03 12.46±0.03 

10-10515 82.55167 1.80581 05301240-1-0148214 3.9±3.9 -1.7±3.9 16.47 1.29 0.27 13.69±0.03 12.79±0.02 12.10±0.02 

11.52±0.06 11.16±0.06 10.87±0.06 10.35±0.06 7.71±0.05 4.38±0.05 3.81±0.07 

10- 10679 82.41315 1.82079 05293913-f0149156 -3.8±3.9 5.0±3.9 16.86 0.38 0.09 13.72±0.02 13.03±0.03 12.53±0.02 

11.77±0.06 11.41±0.06 11.11±0.06 10.63±0.06 8.48±0.05 4.05±0.06 3.29±0.07 

11- 8549 82.99254 1.66050 05315820-1-0139383 4.1±3.8 -1.6±3.8 17.75 0.96 0.21 14.40±0.04 13.72±0.04 13.13±0.04 



Note. — See section [B] 2MASS quality flag is 'A' in all three bands for all sources (see section [2. 3| l. A description of how to construct general PTF catalog 
identifiers is given in table [T] 

''Known CTTS CVSO 35 ijBricefio et al.ll2005l . lioOTh 

''Reported as a candidate WTTS bv lMcGehej ll2006h (SDSS J052700.12+010136.8). 
'^Resembles the 'dipper' objects discussed bv lMorales-Calderon et al. I 1I2OIII '). 
''See section l4T2l for a discussion of zero-point accuracy. 
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Fig. 20. — : Light curves for the candidate CTTS (see section [5]). Continued in figures I^TH^ 
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Fig. 21. — : Light curves for tlie candidate CTTS - cont. 
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Fig. 22. — : Light curves for th e candidate CTTS - cont. 8-10691 (center panel) is a known CTTS 
(jBriceiio et al.ll2005l ). 
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Fig. 23. — : Light curves for the candidate CTTS - cont. 9 -3485 (top panel) was previously reported as a 
candidate weak-line d TTS bvlMcGeheel J2OO6I). 9-6886 (center panel) resembles the 'dipper' 
sources discussed bv iMorales-Calderon et al.l (|201ll) . 
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Fig. 24. — : Light curves for tlie candidate CTTS - cont. 
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Fig. 25. — : Light curves for the candidate CTTS - cont. 
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Fig. 26. — : Color-magnitude diagram as for figure [TTl for the CTTS identified in the PTF Orion data 
matched against the 2MASS catalog (see section [6l). Errors in R represent the standard devi- 
ation in the clipped (f st-99th percentile) standard deviation of the light curves to account for 
variability; errors in Ks are taken to be the same in the absence of other information, added 
in quadrature with the intrinsic measurement errors in Kg. 
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Fig. 27. — : SEDs for the 10 CTTS-classified sources for which Spit zer archive data were ava ilable (nine 
new candidates and the known CTTS, 8-10691 (CVSO SS lBriceiio erallboOSilioOTl )). Symbols 
are as for figure [151 and the gray line again represents an extrapolated photosphere. All appear 
to have clear excesses, as expected for such sources. 



7. Summary 

We have presented some of the initial results from the PTF Orion survey based on the most 
strongly variable light curves, relating to eclipsing binary systems and PMS stars, two of the main 
science objectives of the survey. We have identified 82 eclipsing binary systems, of which 37 are 
detached, and 45 are 'close' (mostly contact binaries). Of the close systems, 11-8774 is among the 
shortest-period W-UMa type systems known, lying in an as-yet poorly sampled region of W-UMa 
parameter space. Furthermore, the distance estimate for this system puts it within the range of 
the Orion association, implying it may be very young: it may therefore be particularly interesting 
for constraining evolutionary models of contact binaries. 

Among the eclipsing binaries, we find a total of nine candidate PMS systems, with implied 
ages of ~ 7-10 Myr if we assume Orion 25 Ori and/or Orion OBla association membership. Their 
positions in the color-magnitude diagram, estimated distances (with the exception of one uncer- 
tain estimate), out-of-eclipse variability, and proper motions (with the exception of one), are all 
consistent with their proposed youth and association membership, making them good candidate 
PMS binaries. Such systems are potentially valuable for constraining evolutionary models of both 
binary systems and individual stars. Furthermore, two of these systems are potentially young con- 
tact binaries. If these are genuinely as young as 7-10 Myr, this could have significant implications 
regarding the timescales of mechanisms by which they could have formed. 

A number of the binary systems are also candidate low-mass systems, apparently with M- 
dwarf primaries; this is also a regime that remains poorly sampled to date. With the exception of 
the two contact systems (which are not expected at such late spectral types), most, if not all, of 
the above mentioned candidate young systems appear to have low-mass primaries based on their 
2MASS colors. We also find an additional two detached systems whose JHKg colors are suggestive 
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of low mass, but with distance estimates more consistent with being background stars; and three 
contact (or near-contact) systems which appear to be unusuahy reddened. 

Some of the systems show somewhat unusual light curves. One in particular, 7-5291, has the 
appearance of a distorted contact- or near-contact binary curve, but has a much longer period than 
is found for contact binaries, and its position on the color- magnitude diagram is consistent with a 
young age at the distance of Orion. One possible explanation is that it is a semi-detached system 
with an inflated PMS primary, perhaps with a disk - although the available Spitzer data do not 
seem to show a clear excess. Assuming it is not a contact binary system gives it a distance estimate 
nearer the 25 Ori/Orion OBla association. Alternatively, the light curve might be exp lained by 
a star being occulted by a w arped disk, in an analogue of A A Tau (jBouvier et al.l 120031 ) and the 
'dipper' systems identified by lMorales-Calderon et al.l ()201ll ). Further data and analysis are needed 
for this object. 

One of the detached systems, 9-2980 appears to exhibit stellar pulsations superimposed on the 
eclipsing binary light curve. The pronounced ellipsoidal variation suggests a fairly closely separated 
system. Further data are needed for coverage of the secondary eclipse, but the pulsations could 
potentially allow for valuable asteroseismological measurements of the system to complement those 
that can be obtained from standard light curve fitting and RV measurements. 

Finally, we have identified 16 candidate CTTS systems, 14 of which are n ew to our knowledge, 
1 of which is aheady known (8-10691, CVSO 35. iBriceho et al.l bood . [2OO3), an d 1 of which ha s 
previously been reported as a candidate WTTS (9-3485, J052700. 12+010136.8; iMcGeheel [2OO6I 1 . 
One of the new CTTS, 9-6886, has the appearance of another 'dipper' source. 

Absolute photometric accuracy in the data set is estimated to be around 3%-5% (» 7% for 
chips 5 and 11). The noise floor in the differential photometry sets in at a precision of around 
4mmag for a single exposure, without any attempt at detrending. The differential photometry 
algorithm developed appears to be relatively robust against systematic trends in all but the worst 
cases, although, owing to known instrument issues at the time of these observations, detrending 
the data may reveal more light curves that could be of interest, particularly on chip 4. 

The data reported here represent an initial cut of only the most strongly variable stars, but 
have yielded a wealth of interesting data. Further followup data on the same field taken with 
PTF at the time of writing will enable confirmation of interesting sources, and help to distinguish 
persistent sources of variability (binaries, planets, etc.) from more transient sources (star spots, 
stellar activity, etc.). A finer search of the PTF Orion data at a lower variability threshold will 
explore a much larger sample of sources. Much remains to be mined in the data. 
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